Abelian deterministic self organized criticality model: Complex 

dynamics of avalanche waves. 

Jozef CernalC 

P. J. Safdrik University in Kosice, Institute of Physics, 
Jesennd 5, O4OOO Kosice, Slovak Republic 

Abstract 

The aim of this study is to investigate a wave dynamics and size scaling of avalanches which 
were created by the mathematical model [J. Cernak, Phys. Rev. E 65, 046141 (2002)]. Numerical 
simulations were carried out on a two dimensional lattice L x L in which two constant thresholds 
= 4 and E^J > Ei were randomly distributed. A density of sites c with the threshold E^.^ and 
threshold Ej.^ are parameters of the model. I have determined autocorrelations of avalanche size 
waves, Hurst exponents, avalanche structures and avalanche size moments for several densities c 
and thresholds Ej.^ . I found correlated avalanche size waves and multifractal scaling of avalanche 
sizes not only for specific conditions, densities c = 0.0, 1.0 and thresholds 8 < E^.^ < 32, in 
which relaxation rules were precisely balanced, but also for more general conditions, densities 
0.0 < c < 1.0 and thresholds 8 < E^.^ < 32, in which relaxation rules were unbalanced. The results 
suggest that the hypothesis of a precise relaxation balance could be a specific case of a more general 
rule. 
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I. INTRODUCTION 

Bak, Tang, and Wiesenfeld (BTW) ^ introduced a concept of self-organized critical- 
ity (SOC) to study dynamical systems with spatial degree of freedom. They proposed a 
simple model with conservative and deterministic relaxation rules to demonstrate the SOC 
phenomenon. Manna ^ designed another conservative SOC model in which stochastic re- 
laxation rules were defined. Striving to find common features of the models and to know 
their basic behaviors stimulated many numerical and theoretical studies during the past two 
decades. 

Based on renormalization group calculations, Pietronero et al. [s^ claimed that both 



deterministic [l| and stochastic models [2| belong to the same universality class, i.e. a 
small modification in the relaxation rules cannot change universality class. Chessa et al. 
4] assumed that finite size scaling (FSS) is common property of both deterministic [U] and 



stochastic [2] models. With FSS the avalanche size, area, lifetime, and perimeter follow 
power laws with cutoffs jj]: 



P(x) = (1) 

where P{x) is the probability density function of x, F is the cutoff function, and and 
are the scaling exponents. The set of scaling exponents (tj., D^) defines the universality 



class jij. A SOC model is Abelian if a final stable configuration (see below) does not depend 
on the relaxation order. The BTW model is Abelian, however the M model is also Abelian 
only if we consider probabilities of many stable configurations. 



Based on numerical simulations and an extended set of exponents, Ben-Hur and Biham jo] 
claimed that the BTW and M models cannot belong to the same universality class. A precise 
numerical analysis of probability density functions P{x) led Liibeck and Usadel Q] to the 

4| 

7|] has been achieved by Tebaldi et al. [8]. They found that 
avalanche size probability density functions P{s) do not display FSS but show a multifractal 
scaling i.e. the avalanche size exponent (Eq. [T]) does not apply to the BTW model. 
Karmakar et al. Q proposed a hypothesis that the presence or absence of a precise relaxation 
balance between the amount released by a site and the total quantity which the same site 
receives when all its neighbors relax at once determines the appropriate universality class. 



same conclusion. Essential progress to understand the discrepancy between theoretical 
and numerical conclusions m 
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Based on the precise relaxation balance hypothesis js] Karmakar and Manna proposed 
a flow chart to classify different SOC models into two universality classes i.e. the BTW and 
Manna universality classes. 



T]sh ow transitions from 
ll| . The models 

m 



The probability density functions of avalanche sizes P{s) in Eq. 
multifractal to FSS scaling for certain densities of disturbing sites | 

are stochastic and non-Abelian with unbalanced relaxation rules jol. In this study, I focus 
on verifying an existence of such transitions for the deterministic and Abelian model [l2[ (In 
the orginal paper 12| the model was incorrectly clasified as non-Abelian) with unbalanced 
relaxation rules. The model 13] displays an anomalous increase of the avalanche size area 
exponents (Eq. [T]) for densities near c = 0.01 and thresholds E^.^ > 16. However, the 
cause of this anomalous behavior is not well understood. I assumed that the transition from 
multifractal to FSS scaling could take place for density c < 0.01 and threshold E^.^ > 16, 
because relaxation rules change character from balanced (c = 0.0) to unbalanced (c > 0.0). 
To characterize avalanche size scaling I investigated avalanche wave dynamics 13|, iM] , Hurst 
exponents [3], avalanche structures js] and avalanche size moments 13|. 

The paper is organized as follows. In Sec. [Ill I repeat a definition of the inhomogeneous 
sand pile model [12] . In Sec. IIIII I determine autocorrelations and fluctuations of avalanche 
size waves, avalanche structures and avalanche size moments for several densities 0.0 < c< 
1.0 and thresholds 8 < E^.^ < 32. Sec. [IV] is devoted to a discussion which is followed by 
conclusions in Sec. |Vl 



II. AN ABELIAN DETERMINISTIC AND CONSERVATIVE SELF ORGANIZED 
CRITICALITY MODEL 

The inhomogeneous sand pile model |l^ is defined on a two dimensional lattice of size 
L X L. Each site i has assigned variables E{i) and Ec{i). The variable E{i) is dynamical 
and it represents a physical quantity such as energy, grain density, and etc. The threshold 
Ec{i) is a static value at site i which is defined only once during initialization of simulations. 
The threshold Ec{i) has two values [l^ : 
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m = { (2) 

A; = 2,3,4,..., 

where c? is a dimension and A; is a natural number. The model has two parameters namely 
the density c = n/L"^ and threshold E^^ where n is a number of sites with the threshold 
Ec{i) = E^^ , remaining L"^ — n sites have the threshold E^ = 4. During initializations of 
simulations, n sites with thresholds -E'c(i) = E^^ were picked out randomly and all remaining 
sites had the threshold E^.^ = 4, thus a set of the thresholds {Ec{i)} represents a quenched 
disorder. A stable configuration is defined by a condition E{i) < -E'c(i) for each site i. Let 
us assume that from a stable configuration we iteratively select i at random and increase 
-E'(i)— E{i) + 1. If an unstable configuration is reached i.e. Ej i) > Ec{i) then a relaxation 



starts. The relaxation rules are conservative and deterministic 



12|: 



E(i)^E(i)-5^AE(e), 



(3) 



^(i + e) ^ ^(i + e) + AE{e) 



(4) 



J2^Eie)=E,{i), (5) 

e 

where e is a set of vectors from the site i to its neighbors. The relaxation rules (Eqs. [3ll5]) 
are repeated until the site i becomes stable. If the neighbors of the site i become unstable 
then avalanche can run on. All unstable sites belong to the avalanche. The relaxations 
given by Eqs. [3}|5]are repeated until a stable configuration is reached, i.e. E{i) < Ec{i) for 
all sites i. Stable and unstable configurations are repeated many times. A total number of 
relaxations during one avalanche is an avalanche size s. 



III. NUMERICAL RESULTS 



Numerical simulations were carried out on two dimensional lattices L x L where the 
linear lattice size L was L = 128, 256, 512, 1024, 2048, and 4096. A density c and thresholds 
E^ = 4, E^^ = 8, 16 and 32 were chosen based on the previous results [12] to cover a 



4 



parameter space in which interesting behaviors were expected (Sec. [T]). For example, near 
the density c = 0.01 au^omalous increase of the avalanche area scaling exponent (Eq. 

d]) has been observed Near densities c = 0.0 and 1.0 local relaxation rules change 

character from balanced (c = 0.0 and 1.0) to unbalanced (c > 0.0 and c < 1.0) jo], [lol |. 
thus a transition from multifractal to finite size scaling could take place in the intervals 
c > 0.0 and c < 1.0. Considering these assumptions, I have selected the density c as 
follows: c = a (surroundings of c = 0.0) or c = 1.0 — a (surroundings of c = 1.0) where 
a = 0, 0.001, 0.002, 0.004, 0.008, 0.01, 0.02, 0.04, and 0.08. In addition, to cover the whole 
interval 0.0 < c < 1.0, I added the sample concentrations c = //lO, where / = 1, 2, . . . , 9. I 
have recorded about 10^ avalanches after initializations of simulations in which an avalanche 
dynamics has to reach the SOC state [1]. To qualify a reproducibility of the results all 
numerical simulations were repeated once for each lattice size L x L, concentration 0.0 < 
c< 1.0 and threshold E^,^ . A comparison of these data sets showed that the results are well 
reproducible. 



A possibility to decompose an avalanche into waves 



15| is a significant advantage of 



computer models. Because avalanche wave dynamics [13|, ll^l can provide valuable initial 
information about the character of an SOC model. An avalanche of size s is decomposed 
into m waves with size Sk, where s = J2T=i^k- A time sequence of avalanche waves Sk is 



used to determine the autocorrelation function 
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where time is t = 1, 2, . . . , and the time averages are taken over 5 x 10^ waves. Autocorre- 
lations C(t, L, c) have been analyzed for the time 1 < t < 1000, lattice sizes 128 < L < 4096, 
selected concentrations 0.0 < c < 1.0 and thresholds = 4, E^^ = 8, 16, and 32 (see 
above). The autocorrelations C{t, L = 4096, c) for the biggest lattice size L = 4096 
are shown in Fig. [1] 1 have observed that for the density c = 1.0, the autocorrelations 
C{t, L = 4096, c = 1.0) agree within experimental error with autocorrelations of the BTW 



model ll3| {C{t, L = 4096, c = 0.0)) . I note that the autocorrelations C{t, L = 4096, c = 0.0) 
are not shown in Fig. [H I have approximated the autocorrelations C{t, L, c = 0, 1.0) by 
a simple function C{t,L,c = 0, 1.0) ~ exp(— at) where a is a decay rate [16]. For more 
general conditions, densities 0.0 < c < 1.0 (Fig. [1]), the autocorrelations C{t, L, c) are more 
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L=4096, E"=8 



(a) 



L=4096, E"=16 



(b) 



L=4096, E"=32 



(C) 





FIG. 1: Autocorrelations C{t, L, c) for the linear lattice size L = 4096, for densities 0.001 < c < 1.0 
and for thresholds = 4, (a) Ej.^ = 8, (b) E^J = 16, and (c) E^J = 32. The autocorrelations 
C{t, L, c = 0) are not shown. 



complex functions than for specific densities c = 0.0 and 1.0. I have found that with increas- 
ing time t, the autocorrelations C{t, L, c) are decreasing functions. An unexpected finding 
is the existence of oscillating components of autocorrelations. For all densities 0.0 < c < 1.0 
(Fig. [D^a)) and the threshold E^.^ = 8 their periods are approximately constant. Ampli- 
tudes of oscillating components decrease with increasing density c and time t. At the given 
time t the autocorrelations C{t, L, c) increase if a density c increases, i.e. if C2 > ci then 
Cit, L, ci) > C(t, L, C2). Near densities c = 0.0 and 1.0, the oscillating parts of C(t, L, c) 
disappear. I have observed more complex behaviors (Fig. [1] (b) and (c)) for thresholds 
E^^ > 16 than for the threshold E^^ = 8 ((Fig. [1] (a)). The oscillating components of 
C(t, L, c) (Fig. [Hb)) have longer periods for densities c < 0.08 and threshold E^^ = 16 
than for the threshold E^^ = 8. However, odd periods were split for densities c > 0.08. The 
same effect take place for the threshold E^^ = 32, but a critical density is higher c = 0.40 
(Fig. [U^c)). After splitting the oscillating components, for thresholds E^.^ = 16 and 32, their 
new periods were approximately equal to the period which has been found for the threshold 
El' = 8. 



Stochastic processes are often characterized by Hurst exponents 17|]. To determine the 
Hurst exponent the fluctuation F{t, L) 1^: 
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F{t,L)=[{Ay{tf),~{Ay{t))lY^' (7) 

is used where y{t) = X]fc=i and Ay{t) = y{k + t) — y{k). If a fluctuation F{t, L) scales 
with the time t as F{t, L — )■ oo) ~ then H is the Hurst exponent. 

1 have determined fluctuations and the corresponding Hurst exponents for the lattice 
size L = 4096, selected concentrations c (see above) and thresholds = 4, E^.^ = 8, 16, 
and 32. For all these parameters fluctuations show two scaling regions 1 < ti < 100 and 
1000 < t2 < 10000. The fluctuations F{t, L), for densities c = 0.000, 0.001, 040 and 0.999 
and threshold = 32, are shown in Fig. [2] to demonstrate the existence of two scaling 
intervals. The Hurst exponents Hi{c) and H2{c) as functions of density c and threshold E[f 
are shown for all parameters (densities 0.0 < c < 1.0 and thresholds 8 < E^^ < 32) in Fig. 
[3l I have observed that the exponents Hi{c) and H2{c) depend on the parameters (density 
c and threshold Eq) in a nontrivial manner (Fig. |3]). For densities 0.0 < c < 1.0 and 
thresholds 8 < < 32 functions Hi{c) are bounded by the interval 0.68 < Hi{c) < 0.81. 
Similarly, the functions -f^2(c) are limited by the interval 0.44 < -f^2(c) < 0.56. The exponents 
Hi{c) = 0.80 and H2{c) = 0.50 are approximately constant for densities c > 0.50. I have 
observed anomalous decreases of functions ifi,2(c) near low densities 0.0 < c < 0.01 (Fig. 
E]). In addition, functions of Hi^2{c) have a decreasing tendency if the second thresholds 
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FIG. 3: Hurst exponents Hi{c)and H2{c) as functions of a densities c and thresholds: E[ = 4, (a) 
E^J = 8, (b) ^^-^ = 16, and (c) Ej.^ = 32 for the lattice size 4096 x 4096. The Hurst exponents 
at given density c, Hi{c) and i^2(c) were approximated using the power laws F{t, L) ~ and 
L) ~ t-^aCc) in the intervals t = 1 - 100 and t = 10^ - 10"^. 



EJ increase. Finally, the functions Hi^2{c) are not symmetric around the density c = 1/2, 
i.e. Hi{0.5 - a) Hi{0.5 + a) and H2{0.5 - a) H2{0.5 + a) for 0.0 < a < 0.5 except 
the specific densities c = 0.0 and c = 1.0, where Hi{0.0) = Hi{1.0) and H2{0.0) = H2{1.0) 
within experimental errors. 

Ben-Hur and Biham [Sjproposed to use avalanche structures to demonstrate a difference 
between BTW l| and M 2] models. An avalanche structure consists of clusters of sites with 
equal numbers of relaxations. The BTW model displays rigorous shell-like structures jsl, ^ 
and the M model displays disordered structures 6| with inner holes 9|. I have analyzed 
several avalanche structures (Fig. Hj) of the inhomogeneous sand pile model [l2| to compare 
them with known structures js, 9|. 

I have observed the avalanche structures which resemble the shell-like structure for den- 
sities c = 0.02 (Fig. Il^a), (d), and (g)) and c = 0.98 (Fig. 1 (c), (f), and (i)). However, 
these structures are not exactly shell-like. A clear visible difference is the existence of holes 
in an avalanche, for example see Figs. Ill^a), (b), and (d). The sizes of these holes vary 
from one site (obviously a site with E^.^ > El) to several sites. Mainly, for low densities 
(c = 0.02) an existence of holes is clearly demonstrated in Figs. Ht^a), (d), and (g) where the 
sites with the threshold E^.^ can absorb and relax more energy than surrounding sites with 
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FIG. 4: Avalanche structures on two dimensional lattice 128 x 128 for parameters = 4, (a) 
EU = 8, c = 0.02, (b) E^j = 8, c = 0.5, (c) E[/ = 8, c = 0.98, (d) E[/ = 16, c = 0.02, (e) 
E[I = 16, c = 0.5, (f) EH = 16, c = 0.98, (g) E[J = 32, c = 0.02, (h) E^j = 32, c = 0.5, and 
(i) Ell = 32, c = 0.98. Lattice sites with the same numbers of relaxations are shown by the same 
color (rainbow pseudo-color coding). 

Ej^. If the sites (Ef ) absorb energy then they are well identified as small holes inside the 
avalanche structure (Figs. Hl^a), (d), and (g)). If the sites (E^^) release more energy than 
their neighbors (with threshold E^) then the sites (-E^^) involve instabilities of many sites 
within a certain distance. These sites must relax to be stable (Fig. 111(d)). At high density 
c = 0.98, the sites with the lower threshold E^ are considered for disturbing sites. A site 
with threshold E^ can receive more energy from neighbors than a critical amount of the site 
(threshold E^.). Then the site {El) must relax more times than neighbors (sites with the 
threshold E^,^) to be stable. The disturbing sites are shown as isolated sites in Figs. 111(c), 
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FIG. 5: Scaling of avalanche size moments {s'^) versus the lattice size 128 < L < 4096. The 
exponents q are q = 1.0, 2.0, 3.0, 3.5 and 3.95 from bottom to top. The model parameters are 
E[ = 4, = 32, (a) c = 0.008, (b) c = 0.40 and (c) c = 0.992. 

(f), and (i). The effects of disturbing sites E^.^ and E^. differs, for the small density c = 0.02, 
the disturbing sites Ej,^ can absorb and relax more energy as their neighbors. However, 
for high density c = 0.98, the disturbing sites E^ can only do more relaxations than their 
neighbors sites. 

I have found new avalanche structures which resemble neither shell-like (the BTW mode) 
nor disordered (the M model) 6|, |9| for density c = 0.50 and thresholds 8 < E^^ < 32 Figs, 
m^b), (e), and (h). Their typical feature is the existence of complex clusters in avalanche 
which more resemble percolation clusters. 

The model (Sec. [Tll) displays shell-like avalanche structures as well as the BTW model 
6|, lot] only for the specific densities c = 0.0 and c = 1.0. 

Sometimes we cannot decompose an avalanche into waves, obviously if we study real 
systems, then avalanche moments |l^ are useful. A property x in FSS obeys scaling given 
by Eq. [1] The moments q are jisl : 



x''P{x,L)dx ~ L'-^*?^ 



(8) 



where crx{q) = (g + 1 — Tx)Dx and Xmax ~ L^"^ . I have determined only avalanche size 
moments (s'') versus the lattice size L which are shown in Fig. [5] for densities c = 0.008, 
c = 0.40 and c = 0.992 and for thresholds E^ = A and El^ = 32. The moments (s^) scale 
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L=4096, E"=8 




FIG. 6: Plots d(Js{q)/dq versus exponent q for the lattice size 4096 x 4096, concentrations 0.0 < 
c < 1.0 and thresholds: = 4, (a) E^^ = 8, (b) E^^ = 16, and (c) E^^ = 32. 

with the lattice size L as well as (s'^) ~ L'^^^'^^ (Eq. [8]), thus a basic requirement is met to 
determine 

Using the functions CTs{q) I determined the plots das{q)/dq versus the exponent q which 
are shown in Fig. [6] for densities 0.0 < c < 1.0 and thresholds 4 < E^^ < 32. I have 
observed that the functions das{q)/dq are increasing if exponents q increase (Fig. [6]) for 
exponents q > 1.0, densities 0.0 < c < 1.0 and thresholds 8 < E^,^ < 32. Surface cuts 
d<7s{q)/dq \q=i.o and das{q)/dq |g=4.o, for exponents q = 1.0 and q = 4.0, as functions of 
density c are shown in Fig. [7] to demonstrate this increasing tendency. I have found that 
d(Ts{q)/dq \g=4.o> das{q)/dq \g=i.o, da^iq) / dq \g=i,o^ const, and da s{q) / dq |, =4.07^ const. 
(Fig. [7]) for densities 0.0 < c < 1.0 and thresholds 8 < E^J < 32. This implies that functions 
das{q)/dq differ from the function das{q)/dq of the BTW model. However, for the specific 
parameters c = 0.0, c = 1.0 and 8 < E^.^ < 32 relaxation rules are precisely balances and 
the functions das{q)/dq versus q {q > 1.0) are identical within experimental errors with the 
function da,{q)/dq of the BTW model Id, Il3|. 



IV. DISCUSSION 



I have found that autocorrelations C 
an autocorrelation of the BTW model 



t, L, c) (Fig. d]) are more complex functions than 



13| . The autocorrelations exhibit oscillating com- 



ponents (Fig. [T]) which periods and amplitudes depend on both densities 0.0 < c < 1.0 and 
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FIG. 7: The plots das{q)/dq \q=i (errors ±0.01) and das{q)/dq |g=4(errors ±0.03) versus density 
c for the hnear lattice size L = 4096 and thresholds E^^ = 4, (a) E^J = 8, (b) Ej.^ = 16, and (c) 
Ej^^ = 32. 

thresholds 8 < E^^ < 32. The oscillating components are probably caused by a periodicity 
in an avalanche wave sequence. I assume that the periodicity could be a consequence of an 
excessive energy storing and release in sites with the thresholds E^^ , apparently when the 
sites have low concentration c < 0.5. In such conditions, relaxations of these sites trigger 
relaxations of surrounding sites, i.e. all surrounding sites within a certain distance from a 



disturbing site have to relax 



12|. This hypothesis could be supported by finding that for 



the same time t amplitudes of oscillating components are decreasing if densities c increase 
and oscillations disappear near the density c = 1.0 (Fig. [1]). The periods were longer for 
low densities c < 0.4 and thresholds E^' > 16 than for the threshold E^^ = 8. This could 
be connected with the stronger effect of disturbing sites E^^ > 16 which can store and re- 
lease much more energy than disturbing sites E^^ = 8. However, a nontrivial dependence of 
periods on thresholds E^^ > 8 and densities c (Fig. [T]) and the cause of period splitting for 
thresholds E^^ > 16 and critical densities c (Fig. [T](b) and (c)) are not understood. For spe- 
cific densities c = 0.0 and c = 1.0, the autocorrelations C(t, L, c = 0.0) and C(t, L, c = 1.0) 
agree well with autocorrelation of the BTW model [9|, 113| . I have found correlated avalanche 
waves under more general conditions in which the relaxation rules are unbalanced ^| for 
densities 0.0 < c < 1.0 and thresholds 8 < E^^ < 32. This is completely opposite to the 
result the hypothesis predicted [9|, HO]. I think that the hypothesis about a precise relaxation 
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balance [9|, Il0| is valid only for the specific densities c = 0.0 and c = 1.0 and thresholds 
8 < E^J < 32. 

Two scahng intervals of fluctuations F(t, L = 4096) (Fig. [3]) support correlated avalanche 



waves 



13| for densities 0.0 < c < 1.0 and thresholds 8 < E^^ < 32. The fluctuations 



F{t, L = 4096) agree well with a fluctuation of the BTW model only for the densities c = 0.0 
and c = 1.0 [sl, [l3| when relaxation rules are precisely balanced js]. For all other parameters, 
densities 0.0 < c < 1.0 and thresholds 8 < E^^ < 32, relaxation rules are unbalanced and 



the hypothesis j9|, [lOj predicts a single scaling region with the Hurst exponent H = 1/2. 



However, I have not found single scaling regions of F(t, L = 4096). So the fluctuations 

nn 

F(t, L = 4096) contradict the hypothesis of precise relaxation balance [9|, [lOj. Asymmetries 
of functions Hi 2(c) to permutations of densities c (Fig. E]) could be a consequence of different 
local effects of disturbing sites E^.^ near density c = 0.0 and disturbing sites E^ near density 
c = 1.0. The Hurst exponents Hi{c) are limited by the interval 0.68 < Hi{c) < 0.8 near the 
density c = 0.0 thus avalanche waves are less correlated than in the BTW model {Hi = 0.8) 
and they are more correlated than in the M model {H = 0.5). The Hursts exponents H2{c) 
are limited by the interval 0.44 < H2{c) < 0.56 which indicate that local perturbation effects 
can change a long-term persistence (antipersistence) [17]. The changes of Hurst exponents 
Hi^2{c) with density c and threshold E^^ demonstrate that local perturbation effects could 
change a global wave dynamics. 

I assume that avalanche wave dynamics (autocorrelations and fluctuations) on a finite size 
lattices L x L, for L > where is a critical length, can provide basic information about 
correlated (the BTW model) or uncorrelated (the M model) nature of waves in avalanches 
not only for the finite size L but also for the size L which goes to infinity. I assume that 
L = 4096 is greater than Lc- Then conclusions regarding an avalanche wave dynamics for 
the finite lattice size L = 4096 could be extended for infinite systems. 

A comparison of avalanche structures, for densities 0.0 < c < 1.0 and thresholds 8 < 
E^^ < 32, shows that the avalanche structures of the model (Sec. HI]) are more disordered 
than shell-like structures of the BTW ra, ^1 model but they are more ordered than structures 
of the M model [6|, |9[. Only for the specific densities c = 0.0 and c = 1.0, the avalanche 
structures are exactly shell-like as well as in the BTW model js, 9|. I assume that connection 
between avalanche structures and autocorrelations near the density c = 0.5 are more weak 
for the model (Sec. [IT]) than for the BTW or M models js, 9|. Because avalanche structures 
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are more disordered (Fig. H]) than structures of the BTW model and more ordered than 
structures of the M model, however avalanche waves are correlated as in Fig. [1] 

A hypothesis about a precise relaxation balance 9| for unbalanced relaxation rules pre- 
dicts FSS where das{q)/dq = const, for all exponents q > 1.0. In such case, an expectation 
d(Ts{q)/dq |g=4= d(Ts{q)/dq \g=i is correct. However, avalanche size moments do not sup- 
port this expectation, because I have found that das{q)/dq |g=4> d(7s{q)/dq \q=i and 
functions das{q)/dq |q=4 and das{q)/dq \q=i (Fig. [7]) depend on a density c for all densities 
0.0 < c < 1.0 and thresholds 8 < < 32. This is a reason why functions dcTs (g) / dq of the 
model (Sec. [Tll) cannot be identical with the function das{q)/dq of the BTW model 

m 

(c = 0.0 and 1.0) for densities 0.0 < c < 1.0. Although, I conclude that avalanche waves 
are correlated in avalanches and avalanche sizes scale as a multifractal for all investigated 
parameters. 

For specific densities c = 0.0 and 1.0, the basic assumptions Hi — H'- = of the precise 
relaxation balance hypothesis {q] are correct. However, all results (Sec. IIIII) show a multi- 
fractal scaling of avalanche sizes not only for the specific conditions but also for the more 
general conditions 0.0 < c < 1.0 and 8 < E^^ < 32 in which \Hi — H'^I = An, < n < 7 
where n is a natural number. Thus I conclude that the hypothesis about precise relaxation 
balance [oj] is valid only for the specific parameters c = 0.0 and 1.0 and E^^ > 4. 

I have not observed a transition from multifractal scaling to FSS of avalanche sizes which 
has been expected for low densities near c = 0.0 (Sec. [11). However, the autocorrelations, 

nnn 

Hurst exponents and avalanche size moments support the previous conclusions pjj, 112|, 116| 
that multifractal scaling is very sensitive to local perturbations for densities near c = 0.0. 

Does the model belong to the BTW class? The answer is not uniform. The model shows 
common features with the BTW model for example, correlated waves in avalanches and 
multifractal scaling of avalanche sizes. However, I have demonstrated that the autocorrela- 
tions (Fig. [1]) are complex functions. Hurst exponents are functions of densities c (Fig. [3]) 

and holes can be found in avalanches while all these features are not typical for the BTW 

□ 

model. A possible difference between the inhomogeneous sand pile model and the BTW 
fl 

model [l[ is supported by functions das{q)/dq (Figs. [6] and [7]) which are not identical with 
the function das{q)/dq of the BTW model except the specific densities c = 0.0 and 1.0. 
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V. CONCLUSION 



Applying the classification scheme 10| on the inhomogeneous sandpile model model (Sec. 
[IT|) places the model in the M universality class. However, I have demonstrated (Sec. IIIip 



that the model belongs neither in the M universality class nor in the BTW universality class 
(Sec. II V|) . I assume that it belongs in a multifractal universality class js] which is more 
general than the BTW class [16]. Based on the wave autocorrelations, fluctuations, avalanche 
structures and avalanche size moments, I conclude that an avalanche wave dynamics and 
avalanche size scaling depend on local relaxation details. In addition, a hypothesis about 
precise relaxation balance jo], lol could be a specific case of a more general rule. The reason 
why a multifractal scaling is very sensitive to local perturbations is not well understood. 
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